syms x;
m=1025;
f=exp(-m*x^2);
vpa(int(f,x,0,inf),5)
